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Abstract. 

The PVLAS collaboration is presently assembling a new apparatus (at the INFN 
section of Ferrara, Italy) to detect vacuum magnetic birefringence (VMB). VMB is 
related to the structure of the QED vacuum and is predicted by the Euler-Heisenberg- 
Weisskopf effective Lagrangian. It can be detected by measuring the ellipticity acquired 
by a linearly polarised light beam propagating through a strong magnetic field. Using 
the very same optical technique it is also possible to search for hypothetical low- 
mass particles interacting with two photons, such as axion-like (ALP) or millicharged 
particles (MCP). Here we report results of a scaled-down test setup and describe 
the new PVLAS apparatus. This latter one is in construction and is based on a 
high-sensitivity ellipsometer with a high-finesse Fabry-Perot cavity (> 4 x 10 5 ) and 
two 0.8 m long 2.5 T rotating permanent dipole magnets. Measurements with the 
test setup have improved by a factor 2 the previous upper bound on the parameter 
A e , which determines the strength of the nonlinear terms in the QED Lagrangian: 
^(pvlas) ^ 2 3 x io~ 21 T~ 2 95% c.l. Furthermore, new laboratory limits have been 
put on the inverse coupling constant of ALPs to two photons and confirmation of 
previous limits on the fractional charge of millicharged particles is given. 
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1. Introduction 

In the absence of matter, Maxwell's equations can be obtained from the classical 
electromagnetic Lagrangian density £ C \ ( m S.I. units) 
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2/x \c 2 'J (1) 
where /i is the magnetic permeability of vacuum and c is the speed of light in vacuum. 
A quadratic Lagrangian leads to linear partial differential equations for the fields, and 
the superposition principle holds, thereby excluding light-light scattering and other non 
linear electromagnetic effects in vacuum. 

With the introduction of Dirac's equation for electrons and Heisenberg's 
Uncertainty Principle, Euler, Heisenberg and Weisskopf in 1936 [1] derived an effective 
Lagrangian density which leads to electromagnetic non linear effects even in vacuum. For 
photon energies well below the electron mass and fields much smaller than their critical 
values, B < B crit = m 2 e c 2 /eh = 4.4 x 10 9 T, E < E ciit = m 2 e c 3 /eh = 1.3 x 10 18 V/m, 
the Euler-Heisenberg- Weisskopf effective Lagrangian can be written as 
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h/mc being the Compton wavelength of the electron, a 



(2) 
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e 2 /(hcA7re ) the 



fine structure constant, m P the electron mass. 





a) Leptonic e*e LbL scattering 



b) Leptonic e + e vacum birefringence 





c) Leptonic e*e vacuum birefringence with second order radiative corrections. 





d) LbL hadronic scattering with 
gluons in the qqbar bubble 



e) Birefringence due to virtual 
spin zero bosons (e.g. axions) 



Figure 1. Feynman diagrams for four field interactions. 



This Euler-Heisenberg- Weisskopf Lagrangian allows four field interactions that can 
be represented, to first order, by the Feynman diagrams shown in Figure [T] a) and 
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b). Figure [T] a) represents light by light scattering whereas Figure [T] b) represents 
the interaction of real photons with a classical field leading to vacuum magnetic 
birefringence. 

Also Quantum Electrodynamics (QED) predicts non-linear effects in vacuum 
leading to birefringence and light-by- light (LbL) scattering through the four-photons 
box diagram [21[3llIJ[5l[6l[71|8j[9]. Furthermore, magnetic birefringence (see Figure 
[He)) and magnetic dichroism could be generated by hypothetical axion-like particles 
( ALPs) pU HI US 13] or millicharged particles (MCPs) pH CGS CEl HO EE] • Finally, the 
coupling of four photons through qq fluctuations is also possible (see Figure[T]d)). In view 
of the values of the masses of even the lightest quarks (much heavier than the electron 
mass) this last contribution should be very small. However, theoretical predictions have 
a very large uncertainty - they span three orders of magnitude [19] ; moreover, it is not 
possible to evaluate the QCD contributions from indirect measurements, unlike what 
happened for the muon g — 2 [201 121] . 

Vacuum magnetic linear birefringence and light-light interaction in vacuum at very 
low energies have yet to be observed. Several experimental efforts are underway to 
detect such effects [221 EHl EH EHl ESI EZl EH EHl EQ] . Before the present work, the 
previous bound on four photon interaction was set by the PVLAS collaboration [31] 
with an upper bound on vacuum magnetic birefringence An^ PVLAS ^ measured with a 
magnetic field B = 2.3 T 

An (PVLAS) <10x 10 -i9 @ 2.3 T and 95% c.l. (4) 

The measurement has been performed at A = 1064 nm, corresponding to a photon 
energy fiw = 1.17 eV. The limit translates in an upper bound for unpolarised photon- 
photon interaction [HI El El El IH] as 

a(7 LAS )(1.17 eV) < 4.6 x 10" 62 m 2 . (5) 

The predicted QED value of the vacuum magnetic birefringence (see below) An' EHW ' 

and photon-photon elastic scattering cross section 

An (EHW) = 2.1 x 10- 23 @ 2.3 T (6) 

a™(1.17 eV) = 1.8 x 10~ 69 m 2 (7) 

1.1. Electrodynamics 

To calculate the magnetic birefringence of vacuum one can proceed by determining the 
electric displacement vector D and the magnetic intensity vector H from the Lagrangian 
density £ehw of Eq. (TSJ) by using the constitutive relations [3] 

D = ^ and H = -^. (8) 
dE dB 
From these one obtains 



D =e E + e A e 
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Let's consider a linearly polarised beam of light propagating perpendicularly to an 
external magnetic field B ext . With D and H, and using Maxwell's equations in media 
one can describe light propagation in an external field. It is evident that the equations 
for light propagation will no longer be linear due to the non linear dependence of D 
and H on E and B, respectively. We substitute E = _E 7 and B = S 7 + B ext in (jHJ) and 
([TO]) , where the index 7 refers to the light, assuming B eXt ^> B 1 . One finds the following 
relations for the relative dielectric constant and magnetic permeability of vacuum: 

4 EHW) =i+io^ e 2 xt r e r w) =i-4A e E e 2 xt 

/T W) = 1 + 4A e 5L ^ = 1 + (11) 

where the indices || and _L refer to light polarisation parallel and perpendicular to B ext , 
respectively. 

From these sets of equations two important consequences are apparent: the velocity 
of light in the presence of an external magnetic field is no longer c and vacuum is 
birefringent with 

An( EHW > = n[, EHW) - nf HW) = SA^lt- (12) 

Numerically this leads to the value given in equation ([6]). In any respect, magnetised 
vacuum behaves like a uniaxial crystal. 



1.2. Post-Maxwellian generalisation 



It is interesting to generalise the non linear electrodynamic Lagrangian density 
correction, which must be parity conserving, gauge and Lorentz invariant, by introducing 
two free parameters 771 and 772: 
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(13) 



where £ = l/_B 2 rit , and rji and 772 are dimensionless parameters depending on the model. 
With such a formulation the birefringence induced by an external transverse magnetic 
field is 

A« (pM) = 2£(V2 - Vi) Bt t . (14) 

This expression reduces to Eq. (Tl2l) if r/i = a/457r and r] 2 = 7/4r/i. It is apparent 
that n\\ depends only on 771 whereas n_|_ depends only on r] 2 . It is also noteworthy that 
if 7]i =772, as is the case in the Born-Infeld model [32], then there is no magnetically 
induced birefringence even though elastic scattering is present 



1.3. Higher order QED corrections 

Figure [1] c) shows Feynman diagrams for the a 3 contribution to the vacuum magnetic 
birefringence. The effective Lagrangian density for this correction has been evaluated 
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by different authors [3S1 ESI EI] and can be expressed as 
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(15) 



This Lagrangian leads to an extra correction An' rad ' to the vacuum magnetic 
birefringence given in equation (fl~2l) 

25a. 



An (rad) 



0.0145 x 3A e B 2 xt 



(16) 



resulting in a 1.45 % correction to An( EHW ). 



l.Jf.. Other hypothetical effects 

Two other important hypothetical effects could cause n^l and An ^ in the presence 
of an external magnetic field transverse to the light propagation direction. They can be 
due either to neutral bosons weakly coupled to two photons, like the axion [TOl [Til [T2| [T3] 
or to millicharged particles (MCP) [H US EH El HE] . 



i.^.i. Axion-like particles (ALP's) Search for axions or axion-like particles using 
laboratory optical techniques was experimentally pioneered by the BFRT collaboration 
[38] and subsequently continued by the PVLAS collaboration with an apparatus installed 
at INFN National Laboratories in Legnaro (LNL) [2211231 EH- In 2006 the detection of 
a dichroism induced by a magnetic field in vacuum was published |39j. Such a result, 
although in contrast with mo del- dependent interpretations of the results of the CAST 
experiment [40J, could have been due to axion-like particles. Soon after, this result was 
excluded by the same PVLAS collaboration [231 EI] after a series of upgrades to the 
apparatus, and almost simultaneously the axion-like interpretation was excluded by two 
other groups [4TJ, [421 S3] with regeneration-type measurements. 

The Lagrangian densities describing the interaction of axion-like particles with two 
photons, for convenience expressed in natural Heaviside-Lorentz units, can be written 
as 

C a = g a( j) a E-B and C s = g s (p s [E 2 — B 2 ) (17) 

where g a and g s are the coupling constants of the pseudoscalar field a or scalar field 
S , respectively. Therefore, in the presence of an external uniform magnetic field B ext 
the component of the electric field of light E y parallel to B ext will interact with the 
pseudoscalar field. For the scalar case the opposite is true: an interaction is only 
possible for the component of E 1 normal to B ext . 

For photon energies above the mass m ayS of such particle candidates, a real 
production can follow: the oscillation of photons into such particles decreases the 
intensity of only one of the polarisation components. On the other hand, even if the 
photon energy is smaller than the particle mass, virtual production can take place, 
causing a phase delay of one component with respect to the other. The phase difference 
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A(p = (j)u — 0j_ accumulated in an optical path L inside the magnetic field region will 
generate a birefringence 

An = |f * (18) 

271 L 

We use the complex index of refraction, written as n = n + in, where n is the 
index of refraction and k is the extinction coefficient, which is related to the absorption 
coefficient /j, by fi — Attk/X. A dichroism can be described by the difference in extinction 
coefficient Ak = k\\ — n± of the medium for the two orthogonal polarisations. In the 
pseudoscalar case, where njj > 1, k?, > 0, rij_ — 1 and «® = 0, and in the scalar case, 
where n s ± > 1, Kj_ > 0, nf, = 1 and kJ = 0, the dichroism Ak and the birefringence An 
can be expressed, in both the scalar and pseudoscalar cases, as PHI EU 155] : 

g„ s B cxt L\ 2 /sinx x 2 



\Ak\ = /cjj = k s ± = 2 f - j l^^J ( 19 ) 

i a i « 1 » 1 ^^ext A sin2a;\ 

|An| = n„ - 1 = n ± - 1 = 1 - — — (20) 



2m l, s V 2x , 

where, in vacuum, x = —^ ± , w is the photon energy and L is the magnetic field 
length. The above expressions are in natural Heaviside-Lorentz units so that 1 T 
^ = 195 eV 2 and 1 m = f = 5.06 x 10 6 eV" 1 . 
In the approximation for which x C 1 (small masses) expression ( 120]) becomes 

A I a i s 1 fcj.s ext' a.s /01 \ 

^ = n?\ — 1 = - 1 = — ! ■ — (21) 

whereas for x ^> 1 

2 ^2 

|An| = nl - 1 = n° - 1 = ext . (22) 

It is interesting to note that the birefringences induced by pseudoscalars and scalars are 
opposite in sign: n| > n" = 1 whereas njj = 1 < n s ± . The different behavior of — 1 
and Um — 1 with respect to L in the two cases where x < 1 and x 3> 1 leaves in principle, 
if the magnet length can be varied, a free experimental handle for distinguishing between 
the two scenarios: this is in fact the case of the new PVLAS setup. 



1.4-2. MCP Consider now the vacuum fluctuations of particles with charge ±ee and 
mass m e as discussed in references [H] and [15]. The photons traversing a uniform 
magnetic field may interact with such fluctuations, resulting in both a pair production 
and a phase delay if the photon energy u > 2m e , and only a phase delay if u < 2m e . 
We are discussing here only the birefringence and hence the (real) index of refraction. 
Moreover, we consider separately the case of fermions and of bosons. 
- Dirac fermions 

Let us first consider the case in which the millicharged particles are Dirac fermions 
(Df). As derived by Tsai and Erber in reference [16], the indices of refraction of photons 
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with polarisation respectively parallel and perpendicular to the external magnetic field 
have two different mass regimes defined by a dimensionless parameter \ (S.I. units): 

3 hio eeB cxt h 



X 



m 2 c 2 



2 m e c 2 

It can be shown that [HJ [T7] 

nZ=l + I?U X )A t B i 
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945^/221/3 (r(|)y 

"7 2 



r(|) 



-4/3 



(3)n,(2) 



for x ^ 1 
for x > 1 



and 



A r 



2 e 4 a 2 A 3 



(23) 



(24) 



(25) 



(26) 



45/x m e c 2 

in analogy to equation ([3]). In the limit of large masses (x «C 1) this expression reduces 
to (fTTl) with the substitution of ee with e and m e with m e in equation (j24p . The 
dependence on 5 ext remains the same as for the well known QED prediction. 

For small masses (x ^> 1) the index of refraction now also depends on the parameter 



X 4 / 3 resulting in a net dependence of n with rather than 5 2 



ext ' 



cxt 



(27) 



cxt 



for x ^ 1 
for x > 1- 



In both mass regimes, a birefringence is induced: 
An Df = ( X ) - i? f (*)] ^ - 

9 45 ^/3 (r (|)) 2 

72 r(i) " 

- Scalar particles 

Very similar expressions to the Dirac fermion case can also be obtained for the 
scalar (sc) charged particle case [IH[T5]. Again there are two mass regimes defined by 
the same parameter x of expression fT23|) . In this case the indices of refraction for the 
two polarisation states with respect to the magnetic field direction are 



n II 



ext 



with 



9 45 7^21/3 (r(|))' 



14 2 r (|) 
The magnetic birefringence is therefore 
An sc = [if (*) - il c (x)] A e B, 

j2 



X 



-4/3 



2 

ext 



(28) 

for x "C 1 
for x > 1- 

(30) 



-A B' 



9 45 tt 1 / 2 2 1 / 3 (r(|))- 

14 2 r (I) 



- 4/ ^e5e 2 xt 



for x ^ 1 
for x > 1- 



(29) 
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As can be seen, there is a sign difference in the birefringence An induced by an external 
magnetic field in the presence of Dirac fermions with respect to the case in which scalar 
particles exist. 



2. Experimental method 

A birefringence An induces an ellipticity ip on a linearly polarised beam of light given 
by 

= ^ S in2$ = 7r^^sin2$ (31) 
2 A 

where L is the optical path length within the birefringent region and A is the wavelength 
of the light traversing it; i? is the angle between the light polarisation and the slow axis, 
defined in our case by the magnetic field direction. The optical path length can be 
increased by using a Fabry-Perot cavity; in this case we think of an effective path length 
L e ff and use a capital letter \1/ to indicate the total induced ellipticity. In fact, given 
a birefringent region of length L within a Fabry-Perot cavity of finesse J 7 , the effective 
path length is L cS = 2J 7 L/tc (see section C below). Today finesses T > 400000 can be 
obtained. 

In the Euler-Heisenberg-Weisskopf case, An depends quadratically on the magnetic 
field B. High magnetic fields can be obtained with superconducting magnets but, as 
we will see below, it is desirable to have a time dependent field either by modulating 
its intensity, thereby changing An or by rotating the field direction, thereby changing 
d. This makes superconducting magnets far less appealing than permanent magnets 
which, today, can reach fields above 2.5 T. In the BFRT experiment the first scheme 
was adopted [38] , whereas PVLAS chose to rotate the magnetic field direction and uses 
now permanent dipole magnets, which are relatively inexpensive to buy, have no running 
costs and have in principle 100% duty cycle, allowing very long integration times. 

We are working with a Nd:YAG laser emitting radiation at 1064 nm. Frequency 
doubled versions exist and could double the induced ellipticity, but at the moment the 
highest finesses have been obtained for infrared light. 

The expected ellipticity must be extracted from the noise within the maximum 
available integration time. For this reason, the magnetic field is made time dependent, 
to move away from dc and limit 1/f noise. Homodyne and heterodyne detections are 
two possible detection schemes. This second technique has been adopted in the PVLAS 
apparatus. 



polariser magnetic field E || i pt , c " ^d u | nto r 

^mh- EH] □ EH ' 

mirror at 2c o Mog mirror n (t) at 0) Mod 



Figure 2. Scheme of the PVLAS ellipsometer. 
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A scheme of the PVLAS ellipsometer is shown in Figure [2J The input polariser 
linearly polarises the laser beam of intensity Iq which then enters the sensitive region 
delimited by the Fabry-Perot cavity mirrors. The laser is phase locked to this cavity. 
After the cavity the laser beam passes through an ellipticity modulator which adds a 
known time dependent ellipticity n(t) to the beam. This modulator ellipticity adds to 
the ellipticity \l/(i) acquired within the magnetic field region. After the modulator, the 
beam passes through the analyser which selects the polarisation perpendicular to the 
input polarisation and finally a photodiode detects In-- The sought for information can 
be extracted from the Fourier spectrum of and from the value of the intensity J out 
before the analyser. 



2.1. Estimate of the effect 

To better understand what follows, it is useful to present some numerical values 
of the different quantities involved in the PVLAS experiment. Considering the 
vacuum magnetic birefringence due to the Euler-Heisenberg-Weisskopf Lagrangian, let 
us determine the ellipticity we expect in the apparatus under construction. The magnets 
have a total magnetic field length L = 1.6 m with a measured field intensity |-B ex t| = 2.5 
T resulting in J B^ xt dL = 10 T 2 m. In the estimate we use the the maximum finesse 
value we have reached: T = 414000. Putting these numbers together leads to 

*pvlas = 2? 3Aj f xtCLL = 3 x 10- 11 . (32) 

Assuming a maximum integration time T max = 10 6 s and a signal to noise ratio SNR = 
1 implies that the sensitivity must be 

spvlas < ^pvlasv 7 ^ = 3 x 10~ 8 Hz" 1 / 2 . (33) 

In principle this is well above the shot noise limit (see formula (|48|) below). 
Unfortunately, as will be briefly discussed below, other noise sources are present which 
limit the ellipsometer sensitivity. We now discuss in detail various issues of the 
measurement technique and of the setup. 



2.2. Heterodyne technique 

Considering the coherence of the light source, a full treatment of the system can be done 
with the Jones matrix formalism |44J. For the purpose of our discussion let the laser 
beam propagate along the Z axis and let the incoming (linear) polarisation define the 
X axis (Figure [3]). The Jones matrix for a uniaxial birefringent element of thickness L 
is given by 

BRF(0)=( 1 + ^ C0S f # sin 20 \ 

v ; I z?/>sin2i? 1-«Vcos2t? J K J 

where ip = Atp/2 (ip -C 1) is now the maximum ellipticity acquired by the light and i? 
represents the angle between the slow axis (n\\ > n±) of the medium and the X axis. 
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Figure 3. Reference frame for the calculations using the Jones matrix formalism. The 
birefringent medium has a thickness L. 

Let E in = -Eu(o) be the entrance electric field; after the magnetic field region one 

has 

E B rf = E ■ BRF • ( ) = Eq\ + 

\0 J \ lip sin 2v 

Assuming no losses and adding the analyser — a polariser crossed with respect to the 
entrance polariser, with Jones matrix A = ( jj ^ j — the intensity will be 

I Tl = I \iip sm2$\ 2 . (35) 

Given the predicted value, the resulting output power, proportional to ip 2 , is 
immeasurably small. 

By introducing an ellipticity modulator (in our case a Photoelastic Modulator, 
PEM), which adds a known sinusoidal ellipticity r](t), the detected signal is linearised 
in ip. In fact, the Jones matrix for the modulator is the same as BRF with $ set at an 
angle of 7r/4 (ip <C rj <C 1): 

MOD = ( * Z f } ) . (36) 

The resulting vector describing the electric field after the modulator will be 

£mod = £u MOD BRF- Q = E ^ ^ _ ^(^ cos^ ) 

Neglecting second order terms, the intensity after the analyser will be 

J Tr (t) = I \ tTJ (t) + sin2^| 2 ~ I [r](t) 2 + 2r)(t)ip sin 2$] (38) 

which now depends linearly on the ellipticity ip. To complete the discussion, one finds 
experimentally that static and slowly varying ellipticities, indicated as a(t), are always 
present in a real apparatus and that two crossed polarisers have an intrinsic extinction 
ratio a 2 , mainly due to imperfections in the calcite crystals. Furthermore, losses in the 
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system reduce the total light reaching the analyser, which we indicate as J out . Therefore, 
taking into account an additional spurious ellipticity term a(t) (since a,^,f)<l these 
terms commute and therefore add up algebraically) and a term proportional to a 2 , the 
total intensity at the output of the analyser will be 

H(t) = /out [cr 2 + \ir]{t) + #sin2tf + %a{t)\] ~ 

^ /out [<J 2 + V(t) 2 + a(t) 2 + 27}{t)1> sin 2$ + 2^(*)a(t)] . (39) 

In PVLAS, to be able to distinguish the large term rj(t)ot(t) from the term 
T](t)if} sin 2$, ■0sin2i? is also modulated in time by rotating the magnetic field direction 
(varying The final expression, explicitly indicating the time dependence of for 
the intensity at the output of the analyser is therefore 

/■&(*) = /out [c^ 2 + r/W 2 + a(t) 2 + 2r){t)ij;sm2#{t) + 2 V (t)a(t)] . (40) 



2. 3. Optical path multiplier 

To increase the ellipticity induced by the birefringent region of thickness L one can 
increase the number of passes through it. Either a multi-pass cavity or a Fabry-Perot 
cavity can be used for this purpose. In the PVLAS experiment a Fabry-Perot has been 
chosen. In a multi-pass cavity the induced ellipticity is multiplied by the number of 
passes iVp ass through the region. With a Fabry-Perot cavity the analogy to a multi-pass 
cavity is not immediate since one is dealing with a standing wave. 

Let t, r be the transmission and reflection coefficients, and p the losses of the mirrors 
of the cavity such that t 2 + r 2 + p 2 = 1. Let d be the length of the cavity and 5 = 4ird/X 
the roundtrip phase for a beam of wavelength A. Then the Jones matrix for the elements 
of the ellipsometer after the entrance polariser is 

oo 

ELL = A SP MOD ■ t 2 e lS/2 [ERFVe^] B ■ BRF (41) 

n=0 

where SP describes spurious ellipticities. Because r 2 < 1, ELL can be rewritten as 

ELL = A SP MOD • t 2 e tS/2 [I - BRF 2 rV 5 ] ~* ■ BRF (42) 

with I the identity matrix. With the laser phase locked to the cavity so that 5 = 27rm, 
where m is an integer number, the electric field at the output of the system will be 



En = £ -ELL-( ^ = E ^- — I I ,■- 1 (43) 

V °t 2 +p \ ia{t)+iri{t)+i T -—j>il)sm2# 1 





1 - 



and the intensity, including mirror losses and polariser extinction, 



/ Tr (t) 



'out 



a 2 + 



1 + r 

ia{t) + irj(t) + % ( ) ip sin 2d 



(44) 



This expression is at the basis of the ellipsometer in the PVLAS apparatus 
Small ellipticities add up algebraically and the Fabry-Perot multiplies the single pass 
ellipticity ?/>sin2i?, generated within the cavity, by a factor (1 + r 2 )/(l — r 2 ) w 2J r /7r, 
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where T is the finesse of the cavity. The ellipticity signal to be detected is therefore 
= (2J7tt)^ sin 2$. 

The Fabry-Perot as an optical path amplifier works well until the accumulated 
phase difference due to birefringences within the cavity is less than tt/2. Beyond this 
value the interferometer has two separate resonances for two orthogonal polarisations 
[15] and therefore elliptically polarised light cannot be resonant. As total birefringence 
also includes the intrinsic birefringence of the cavity mirrors, two mirrors with similar 
characteristics must be used and their birefringence axes aligned. 

2-4- Fourier components 

In the PVLAS experiment, T)(t) = rj cos(27rz/Mod^+6 l Mod) and the magnetic field direction 
is rotated at a frequency % ag : = 27n/ Mag t + '# Mag . A Fourier analysis of the intensity 
lTr(t) of equation (jHJ) results in four main frequency components each with a definite 
amplitude and phase. These are reported in table [U 

Table 1. Intensity of the frequency components of the signal after the analyser. 



Frequency 


Fourier component 


Intensity/ J out 


Phase 


dc 

^Mod 
^Mod ± 2i/ Mag 
2^Mod 


-^dc 

I'Mod 

-^VMod±2t/Mag 


^ 2 + «dc + % 2 /2 

2«dc^o 

IT 
m — y> 

ri/2 


#Mod 
#Mod i 2^Mag 
2#Mod 



The presence of a component at UMod ± 2z/M ag in the signal identifies an ellipticity 
induced by the rotating magnetic field. Furthermore the phase of this component must 
satisfy the value in table [U In the presence of a signal above background with the 
correct Fourier phase, the ellipticity \& = (2J r /7r)^ can be calculated from J out , from 
the Fourier components ^ M od±2v Mag , and from hu Mod as the average of the two sideband 
signals: 

^ ~ | "^Mod+ 2iy Mag _j_ -^Mod - 2 ^Mag j (45) 

/2I out I 2uMod y/2I out l2v Mod J 

2.5. Noise considerations 

Indicating with -R^ Mod ±2^ Mag the noise spectral density at the signal frequencies, and 
assuming Ru Mod +2u m „ e = Ru Mod -2u Mag ^ the sensitivity spectral density s of the ellipsometer 
for a unity signal to noise ratio is 

s _ ^Mod+gi^Mag /^g\ 
^/4I out I 2uMod 
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In principle, the noise limit for such a system is determined by the rms shot noise 
4hot °f the current i& c \ 



^shot 



y/2ei dc Av = J 2eI out q + y + "do) ^ (47) 

where q is the efficiency of the photodetector, and Av the bandwidth. In the case 
t]q ^> a 2 , ot\ c , the dc current will depend only on r] and by substituting -Rj, Mod +2i/ M a g = 
ishot/ilV Af) into Eq. fj46l) . the shot-noise sensitivity spectral density s sno t becomes 



Sshot = J- T • (48) 

Assuming an intensity I ont = 5 mW and the efficiency of the diode q = 0.7 A/W, 
equation above gives s sno t ~ 7 x 10~ 9 Hz -1 / 2 . 

Other intrinsic noise sources are photodiode dark current noise idark which generates 
an output noise signal density Vdark = Gid&rk, with G the transimpedance gain 
of the photodetection system, Johnson rms current noise spectral density in the 
transimpedance resistor in the amplifier of the photodiode ij = ^/Ak B T/G, and 
relative laser intensity current noise spectral density 2rjn = lout Q 4 RIN(V), where 
RIN(z/) = I n oise(^)/% is the relative amplitude noise spectral density of the laser light. 
These noises must be kept below i s hot at a frequency near i/ Mod in order to reach 
the theoretical sensitivity. Complete expressions for these noise contributions to the 
ellipticity spectral noise density can be obtained from Eq. (j4"6|) as functions of the 
modulation amplitude rjo: 

s - h " = Vw(^!^j (49) 

Kiark 1 /rn \ 
Sdark = —pT 7 (50) 

G /out q Vo 
[MsT 1 

Sj = \~F^1 ( 51 ) 

V G i out q rjo 

SMN = RIN( ^ d ) «gg , (52) 

Vo 

3. Experimental studies 

3.1. PVLAS - LNL: main features and limitations 

The previous PVLAS apparatus was set up at LNL in Legnaro, Italy and took data 
from 2002 to 2007. It featured a vertical assembly of the ellipsometer with a 6 m 
long Fabry-Perot cavity, the injection and detection optics installed on different optical 
benches resting on the ground and on the top of a granite tower nearly 8 m high, 
respectively. A superconductive dipole magnet about 1 m long, installed in a cryostat 
resting on a turntable, provided a rotating magnetic field orthogonal to the light 
path in the Fabry-Perot cavity. The magnet was operated with magnetic fields up 
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to 5.5 T. The magnet support and the optics tower had different foundations, and 
were hence mechanically decoupled. Although this apparatus has set best limits on 
magnetic vacuum birefringence and photon-photon elastic scattering at low energies 
[3Tj . it suffered from several limitations: 

• high stray field when operating the superconducting magnet above 2.3 T, due to 
saturation of the iron return yoke; 

• limited running time due to liquid Helium consumption; 

• high running costs for liquid Helium; 

• the granite tower could not be seismically insulated, due to its size and 
configuration; 

• with a single magnet, a zero ellipticity measurement was only possible with field 
off, hence not in the same experimental condition as in the magnetic birefringence 
measurements. 

To solve the first three problems, we have now chosen to work with permanent 
magnets instead of superconducting ones. In permanent magnets the stray field can be 
made much smaller and the duty cycle can be as high as 100%. Seismic noise also limited 
the sensitivity of the LNL setup; instead, a large improvement of this parameter can 
be obtained by installing the whole ellipsometer on a single seismically isolated optical 
bench, as shown in ref. [16]. The last point in the list above is extensively addressed in 
a separate section. 

A completely new setup has been designed according to these lines and is being 
installed in a clean room of the Ferrara section of INFN. The new set up is built 
using a single 4.8 m long granite optical bench actively isolated, two identical rotating 
permanent dipole magnets in Halbach configuration, each 0.8 m in length with magnetic 
field strength B = 2.5 T and 2 cm bore. The magnets design features a magnetic field 
shield strongly reducing stray fields. We plan to rotate them at z/ Mag = 10 Hz, thanks 
to their reduced mass and good balancing. The optical setup features a 2 W tuneable 
Nd:YAG laser frequency- locked to a 3.2 m long Fabry- Perot cavity using a modified 
Pound-Drever-Hall scheme [17] and an entirely non magnetic optomechanical setup. 
The motion of the in-vacuum optical elements is based on stepper piezo motors which 
maintain their position when powered off. The whole ellipsometer is kept in a vacuum 
chamber pumped with non evaporable getters, which are intrinsically non magnetic. 
Two 1 m long glass tubes traverse the bore of the magnets; each tube is pumped at 
both ends. 

3.2. Two magnet configuration 

As reported in ref. [23J, the rotating stray magnetic field and the vibrations associated 
to the masses in rotation can act on the optics and fake ellipticity signals. It is therefore 
desirable to perform measurements of the zero of the ellipticity scale in conditions as 
close as possible to those in which an authentic signal is expected. This is not possible 
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with one single superconducting magnet, because when there is no field in the magnet 
bore, also the stray field is absent. By using two identical magnets and orienting their 
fields at 90° the net ellipticity generated by the magnetic birefringence is zero. Running 
the system with the magnets parallel and perpendicular will allow the identification of 
a real physical signal with respect to some spurious signal due to stray field. 




Figure 4. Optical scheme (upper panel) and photograph (lower panel) of the test 
apparatus in Ferrara. At the centre one can see the two dipole permanent magnets. 
The optics is supported by antivibration stages whereas the magnet supports are on 
the floor. 

This idea has been verified with the scaled down test setup using the Cotton- 
Mouton effect of oxygen jH]. The principle of the test ellipsometer is identical to that 
of the final system in construction. Fig. shows the scheme of the test apparatus and 
Fig. |3b a picture of it. At the centre one can see the two permanent magnets each 
generating a 20 cm long magnetic field of maximum intensity B = 2.3 T. For these 
magnets the measured J B 2 dL is 1.85 T 2 m. The whole optical setup is placed on a 
seismically isolated optical bench, whereas the magnets are supported by a separate 
structure resting on the floor and thereby mechanically isolated from the optics. The 
finesse of the cavity was T = 240000 and the oxygen pressure inside the apparatus was 
0.278 mbar. 
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l„ . 7.5 mW 
V 2 „= 28 mV 
G = 0.7 MV/W 

Oxygen pressure = 0.278 mbar 
Finesse - 240000 
B Z L= 1.85 T 2 m 

An u - 2.06. 10"' Z (T Z m)"' 
Ref [47]: (2.29±0.08)- 1 0~ ,z (T z m)" 




Magnet rotation frequency = 1 .5 Hz 

Parallel magnetic fields 

Perpendicular magnetic fields 



-10 10 

Frequency from carrier signal [Hz] 



Figure 5. Fourier spectrum around the carrier frequency z^Mod- Ellipticity 
measurements with the magnets parallel (red) and with the magnets perpendicular 
(black). The large sidebands at 2^Ma g are due to the Cotton Mouton effect in oxygen 
gas. No error is given on the measured value of An u (for a definition see Ref. [49]), as 
we could not control gas composition. 



Figure [5] shows the ellipticity Fourier spectrum around the carrier frequency 
^Mod- Clear sidebands can be seen at twice the rotation frequency of the magnets 
(^Mag = 1-5 Hz). Shown in red is the spectrum with the magnets in a parallel 
configuration whereas the black curve refers to the magnets perpendicular. The signal 
attenuation factor in the perpendicular configuration with respect to the parallel one 
is about 80: we can conclude that the parameter B 2 eY± L for the two magnets is equal 
to within about 1-2% and, most important, that the principle is correct. The presence 
of the unwanted sidebands at i^Mag indicates that a small spurious component could 
also contribute to the 2z/M ag peak, thus limiting the cancellation of the signal in the 
perpendicular configuration. 

We believe the improvement of using two magnets instead of one will be crucial in 
understanding the ellipsometer, and that such a measurement with two rotating magnets 
whose relative field directions can be changed is imperative. This is a straightforward 
way of guaranteeing the authenticity of a signal. 

3.3. Vacuum ellipticity measurements with the test setup 

With the test setup (Figure H]), featuring two small magnets, measurements have been 
performed to understand its limits and optimise the new apparatus in construction. 
Two different noise sources exist and are under study: wide-band noise and spurious 
signals at the magnet rotation frequency and its harmonics. Below we briefly report 
some results in the two cases. 



Vacuum magnetic birefringence: PVLAS 



17 



3.3.1. Sensitivity - wide-band noise Noise measurements were first performed without 
the Fabry-Perot cavity. In this configuration we successfully excluded noise sources 
from readout electronics and optical elements other than the cavity mirrors, practically 
reaching the expected sensitivity of s no cavity = 6 x 10~ 9 Hz _1//2 , dominated by srin, and 
a noise floor of ^Hoor = 1-5 x 10~ 10 with 1600 s integration time. 




Figure 6. Calculated noise contributions for the PVLAS test apparatus (solid curves). 
The green curve is the Johnson noise contribution from the transimpedance resistor in 
the current to voltage converter, the red curve is the dark current noise contribution 
from the photodiode, the blue curve is the shot noise contribution and the brown curve 
is the relative intensity noise contribution measured at the output of the cavity. The 
circle represents the measured ellipticity sensitivity of the apparatus. 

With the introduction of a cavity with finesse J 7 = 240000 the noise increased to 
■^cavity = 3 x 10~ 7 Hz -1 / 2 at about 6 Hz. This was significantly more than what was 
expected from the reduction of I out due to cavity losses (see Eq. ()48l)). The actual 
ellipticity noise budget situation is depicted in Figure |6] as a function of the modulation 
amplitude rj |31] ; in the figure, Relative Intensity Noise (RIN), Johnson noise in the 
amplifier transimpedance, intrinsic current noise in the photodiode and shot noise are 
summed quadratically. The circle marks the observed noise level. This unexplained 
noise is under study. We suspect variations of the intrinsic birefringence of the mirrors. 

We also observed that the magnet rotation did not contribute to the wide-band 
noise, indicating a good isolation between the magnet support and the optical setup. 

3.3.2. Spurious peaks With the magnets in rotation we often observe ellipticity peaks 
varying from a few 10~ 8 to a few 10~ 7 whereas sometimes such peaks are not present. 
The frequencies of these peaks are at harmonics of the magnet rotation frequency. The 
variability of these peaks from one run to another seems to depend (in a non reproducible 
way) on the adjustment of the input and output polarisers which is done with motorised 
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stages. To study the dependence of such peaks on the magnet orientation, a field probe 
is present near the output side magnet. Changing the relative orientation of the two 
magnets does not change the amplitude of these peaks but does change their phase; the 
result of this study indicates that the more sensitive part of the apparatus is the entrance 
optics. All the motorised stages have small electric motors which might couple to the 
rotating stray magnetic field and may introduce beam jitter and therefore ellipticity. 
The substitution of all these stages is under way. 




Figure 7. Top: ellipticity noise spectrum in a 0.78 Hz frequency band around 2i/Mag- 
The spectrum is demodulated with a lock-in amplifier with respect to the i^Mod carrier. 
The integration time was T = 8192 s. Bottom: histogram of the ellipticity noise values 
plotted above. A fit with the Rayleigh distribution is superimposed with an ellipticity 
standard deviation a = 3.35 x 10 -9 . The arrow indicates the value of the ellipticity 
Fourier spectrum at exactly 2^Mag = 6.5 Hz. 
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3.3.3. Noise floor measurements With the test apparatus in a condition in which the 
peak at 2z/M ag was not present, measurements of a few hours have been done. The 
magnet rotation frequency was z/M ag = 3.25 Hz. In the top panel of Figure [7] we report 
the amplitude Fourier noise spectrum of the ellipticity in a frequency band of 0.78 Hz 
centred at 2z/M ag for the best measurement; phase information is not reported. The total 
integration time was 8192 s acquired in a single time record. 

In the bottom panel of the figure, a 50 bins histogram of the same data is given. A 
vertical arrow in the same figure indicates the ellipticity value in the Fourier spectrum 
bin at exactly 2z/M ag — 6.5 Hz. The probability density function for a noise signal 
with equal standard deviations o for the 'in phase' and quadrature components is the 
Rayleigh function: P(r) = N-^e ^ . A fit with the Rayleigh function is superimposed 
to the histogram. The resulting ellipticity standard deviation is a = 3.35 x 10~ 9 which 
translates, at 95% c.L, in a birefringence limit induced by the magnetic field of 

An < v/-21n(l - 0.95) = 4.6 x 10~ 20 . (53) 

The limits obtained for the parameter A e and the photon-photon total elastic scattering 
cross section (at 95 % c.l.) are: 




0.0001 0.001 0.01 

Boson mass m a s [eV] 

Figure 8. PVLAS existence bounds for ALP's. The excluded region is above the 
curves. For m a < 1.5 meV the best limit is still given by the 2008 rotation measurement 
|23j . Above this mass value a new limit is set by the measurement presented in this 
work. For the sake of comparison, the outdated 2008 ellipticity cutves is drawn (dashed 
red) . 
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MCP boson mass [eV] 



Figure 9. PVLAS existence bounds for MCP's. The excluded regions are above the 
curves. Top: Dirac fermions; bottom: scalar bosons. In both cases the best limits are 
still set by the 2008 rotation measurement. The continuous blue rotation curves are 
taken from Ref. [50]. The four ellipticity curves, drawn according to Eqs. (f27| and 
([30| , correspond to the outdated 2008 limit (dashed red) and to the new ellipticity limit 
(continuous red) . The two branches of each of the ellipticity curves are not connected 
in the mass range around \ = 1 (dotted black line) , where the difference of the indices 
of refraction changes sign. 

As shown in Figure El the ellipticity measurement presented here also sets new 
laboratory limits on the existence of axion-like particles for masses above 1.5 meV. Below 
1.5 meV the best laboratory limit is still given by the PVLAS - LNL 2008 dichroism 
measurement [23]. In the case of millicharged particles we confirm the best exclusion 
limits set by the 2008 rotation measurements (see Figs. [9]) for both fermions and scalars. 
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4. Other experimental efforts 

Three other experiments aimed at studying vacuum magnetic birefringence exist: the Q 
& A experiment, the BMV experiment and the OSQAR experiment. All of them study 
the propagation of a laser beam through a magnetic field region and use a Fabry-Perot 
interferometer to increase the optical path length. A number of other proposals have 
been put forward to measure vacuum magnetic birefringence [521 [26j [51] or other related 
quantities (53] with different techniques: they will not be described here. As far as the 
OSQAR experiment [29] is concerned, the information available in the literature at the 
time this article was written was rather incomplete. 

As shown before, vacuum magnetic birefringence appears in a number of different 
theories. In the case of the Euler-Heisenberg-Weisskopf Lagrangian, An is proportional 
to _Bg Xt through the parameter A e = An/(3B^ xk ). On the other hand, the optical 
quantity the experiments are measuring is the ellipticity generated by the induced 
birefringence. The relation between birefringence and ellipticity is given in formula 
(!3T|) . The three experiments have therefore two distinct parts: the optical ellipsometer 
and the magnetic field source. Figures of merit for the optical setup are the ellipticity 
sensitivity s at the frequency of interest, the signal amplification factor defined by the 
Fabry-Perot cavity finesse and the wavelength used. 



Table 2. Main parameters of the three ongoing experiments to measure vacuum 
magnetic birefringence. Effect modulation frequency / mo d is 2i/Mag for PVLAS and 
Q & A, and 1/T pulsc for BMV. 



Experiment 


PVLAS 


Q & A 


BMV 


Status 


Achieved /Planned 


Achieved/Planned 


Achieved /Planned 


Wavelength (nm) 


1064 


1064/532 


1064 


Magnetic dipole 


permanent 


permanent 


pulsed 


jB 2 dL (T 2 m) 


1.85/10 


3.2/19 


25/600 


Average B cxt (T) 


2.15/2.5 


2.3/2.3 


14/30 


Finesse 


2.4 x 10 5 /> 4 x 10 5 


3 x 10 4 /1 x 10 5 


5 x 10 5 /1 x 10 6 


QED ellipticity (Eq. [3TJ 


3 x 10~ 12 /3 x 10 -11 


7 x 10" 13 /3 x lO" 11 


9 x 10- u /5 x 10" 9 


Detection scheme 


heterodyne 


heterodyne 


homodync 


Effect mod. freq. / mo( i 


6 Hz/20 Hz 


26 Hz 


500 Hz 


Duty cicle D t 


~ 1 


~ 1 


3 x 10" 6 


s & / mod (Hs-V*) 


3 x 10 -7 /3 x 10~ 8 


1 x 10~ 6 /1 x 10 -8 


5 x 10~ 8 /7 x 10~ 9 


Seff © /mod (HZ-V2) 


3 x 10~ 7 /3 x 10" 8 


1 x 10~ 6 /1 x 10~ 8 


3 x 10~ 5 /4 x 10~ 6 


An e ff sensitivity (Hz -1 / 2 ) 


1.7 x HT 18 /2.5 x 10~ 20 


3.0 x 10~ 17 /7.4 x 10~ 21 


2.6 x 10" 16 /3.0 x 10~ 18 


A c sensitivity (T~ 2 Hz -1 / 2 ) 


1.2 x 10~ 19 /1.3 x 10" 21 


1.8 x 10~ 18 /4.7 x 10~ 22 


4.4 x 10" 19 /1.1 x 10~ 21 


Time for SNR = 1 


260 yr/12 d 


63 x 10 3 yr/1.4 d 


3.6 x 10 3 yr/8.3 d 



In Table H] the main features of the ongoing experiments are summarised with an 
estimated running time necessary to reach this goal. As the scheme of the Q & A 
experiment [28] closely resembles PVLAS's one, the interpretation of its parameters 
is straightforward. The BMV experiment deserves a few more words: it employs 
an intense pulsed magnetic field lasting T pmse = 2 ms FWHM (total duration about 
4 ms), with a repetition rate of five shots per hour to allow the magnets to cool 
down. Therefore on average the effective duty cycle with the magnetic field ON is 
D t = 2 ms x (5/3600 s) ~ 3 x 10~ 6 . The analysis correlates the measured ellipticity 
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during each pulse, assumed to be proportional to -B 2 xt , with the magnetic pulse shape, 
also taking into account the dynamical response of the Fabry- Perot whose decay time is 
of the order of the magnetic pulse rise time [51]. The analysis gives Ans = An/£?g ltered 
for each pulse, where -Entered is the effective field filtered by the first order response of 
the Fabry-Perot cavity The BMV collaboration therefore presents results in terms 
of the birefringence per tesla squared per pulse. At present BMV has obtained a 
best Aub = 5 x 10~ 20 T~ 2 per pulse and claims that a future upgrade will able 
to reach SNR = 1 in the QED measurement with 1000 pulses. Given a number of 
pulses iVp U i ses with equivalent noise, the limit on An B will scale with 1/ ^N pu \ ses . In 
the detection of a birefringence much smaller than Ane many pulses will be needed 
and therefore, to compare their method with the Q & A or PVLAS capability, it is 
reasonable to rescale Ann per pulse to the average sensitivity in An achievable in 1 
s: An e ff = Ane -£> 2 xt a/ T pu i sc / D t . Since BMV has a small duty cycle D t their effective 



optical sensitivity reduces to s e ff = s/ \fD~t- 



It is apparent that all the experiments are still far from their goal, the chance of their 
eventual success relying on improvements of the optical sensitivity that, unlike hardware 
upgrades, cannot be given for granted, as they would require complete understanding 
of the noise sources limiting the sensitivity. Indeed the shot noise limited sensitivity of 
all the experiments, given by formula ( HHjl (due to the homodyne technique adopted by 
BMV they have an extra factor of l/y/2), is about a factor between 10 to 100 below 
the current achieved sensitivities. No clear explanation to this discrepancy has yet 
been given. Note that the best value of sensitivity in ellipticity in ellipsometers using 
high finesse Fabry-Perot interferometers has been reported by PVLAS in a small test 
ellipsometer without the implementation of a mag netic field [46] of 2 x 10~ 8 Hz _1/2 above 
5 Hz. The best reported sensitivity of a complete apparatus is by BMV with a sensitivity 
of 5 x 10~ 8 Hz -1 / 2 over a time span of 2 ms. If in Table H] the QED ellipticity planned 
by each experiment is combined with the presently achieved sensitivities, the last line in 
the table would become (3.2, 39, 1.5) yr for the three experiments, respectively, instead 



5. Conclusions 

We have presented the physics the PVLAS experiment is aiming at studying and have 
briefly discussed the experimental method. Noise sources are being studied on a bench- 
top small test apparatus in Ferrara, Italy, in view of the construction of the final 
apparatus with which we hope to measure for the first time the magnetic birefringence 
of vacuum due to vacuum fluctuations. 

We have discussed the importance of using not one but two dipole magnets whose 
relative directions can be made perpendicular to each other in order to have a zero 
effect condition with the magnetic field present. This is vital to study and eliminate 
spurious signals generated by the rotating field. A proof of principle measurement using 
the Cotton-Mouton effect in oxygen gas has been done with the test apparatus. 




of (12, 1.4, 8.3) d. 
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Noise floor measurement results have been presented. A new limit on the parameter 
A e describing non linear electrodynamic effects in vacuum has been obtained: A e < 
3.3 x 10~ 21 T~ 2 . This value, obtained by using two small compact permanent magnets, 
improves the previous limit by a factor 2. 

Finally, the perspectives of the three different experiments based on similar 
techniques active on this subject have been discussed. Ellipticity noise sources are 
under study by all groups so as to improve their optical sensitivity. Without such 
improvements the measurement of vacuum magnetic birefringence is still out of reach. 
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